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A  statistical  mechanical  model  for  adsorption  is  used  to  produce 
voltammograms  with  sharp  spikes  similar  to  those  seen  in  recent  studies  of 
the  underpotential  deposition  of  metals  onto  single  crystal  electrodes.  The 
voltammogram  of  the  underpotential  deposition  of  copper  on  the  (111) 
surface  of  gold,  which  contains  two  spikes,  can  be  reproduced  by  the 
model. 


INTRODUCTION 

Several  recent  studies  of  the  underpotential  deposition  of  metals  onto  single  crystal 
electrodes  have  produced  voltammograms  (1)  containing  sharp  spikes  (2)-(7).  These 
spikes  appear  to  be  associated  with  first  order  phase  transitions  which  occur  on  the  nearly 
perfect  lattice  of  adsorption  sites  on  the  electrode  surface. 

We  have  recently  studied  a  statistical  mechanical  model  for  adsorption  (8)-(12)  in  which 
a  solution  of  ions  is  in  contact  with  a  hard  wall  containing  a  lattice  of  sticky  sites.  The 
model  is  equivalent  to  a  two-dimensional  lattice  gas  in  which  the  adsorbed  atoms  have  an 
activity  which  is  a  product  of  a  stickiness  parameter  and  the  contact  density  of  the  ions  at 
the  electrode  surface.  Lateral  interactions  are  also  present  between  neighboring  adsorbed 
atoms. 

Using  eight  exact  coefficients  from  the  known  high  and  low  activity  series  for  the 
fraction  of  occupied  sites  (13),  (14),  we  have  constructed  approximate  ad.sorption 
isotherms  for  the  model  (1 1),  (12)  which  are  a  natural  extension  of  the  Langmuir  isotherm. 
The  time  derivative  of  the  fraction  of  occupied  sites  gives  a  measure  of  the  intensity  of  the 
current  of  the  voltammogram  accompanying  the  underpotential  deposition  process.  The 
model  voltammogram  contains  sharp  spikes  which  arc  associated  with  first-order  phase 
transitions  which  occur  in  the  lattice  gas  (1 1),  (12). 

In  the  next  section  of  the  present  paper  we  discuss  the  adsorption  model  and  the 
voltammograms  which  are  derived  from  it.  In  the  final  section  of  the  paper  we  use  the 
model  to  reproduce  the  voltammogram  of  the  underpotential  deposition  of  copper  onto  the 
(111)  surface  of  gold.  The  two  sharp  spikes  which  are  observed  in  the  experimental 
voltammogram  (3)  can  be  reproduced  by  our  model  (12). 


THE  MODEL 


We  have  recently  considered  a  statistical  mechanical  model  for  adsorption  in  which  a 
solution  of  ions  is  in  contact  with  an  electrode  surface  which  is  modeled  as  a  smooth,  hard 
wall  with  a  lattice  of  sticky  adsorption  sites  (11).  The  three-dimensional  model  can  be 
shown  to  be  equivalent  to  a  two-dimensional  lattice  gas  of  the  adsorbed  atoms  or  ions.  In 
the  equivalent  lattice  gas,  the  adsorbate  species  has  an  activity  z  which  is  a  product  of  a 
stickiness  parameter  X  associated  with  the  lattice  site  times  the  density  p®(0,  *f0  of  the 
ions  at  the  electrode  surface,  which  has  a  dimensionless  potential  bias 


(11 


Here  p  =  (kT)~^,  e  is  the  electronic  charge,  V'(O)  is  the  potential  at  the  electrode 
surface,  and  V/_  .  is  the  potential  of  zero  charge.  We  approximate  the  dependence  of 
pO(0,  ¥0  on  as 


pO(0,  ¥0=pO(0,0)e-5'^  , 


[2] 


where  ^  is  the  electrovalence  of  the  ion  in  solution. 

The  atoms  of  the  lattice  gas  have  a  first-neighbor  interaction  energy  w(r)  which  is  the 
pair-potential  of  mean  force  of  the  adsorbed  species  at  a  first  neighbor  separation  r  on  the 
electrode  surface.  That  is, 


w(r)  =  -kTlng2(r)  , 


[31 


where  g2{r)  is  the  pair  correlation  function  of  the  adsorbed  species  at  a  distance  r. 

For  the  lattice  gas  on  .several  different  lattices  with  lateral  interactions  extending  only  to 
first  neighbors,  many  coefficients  have  been  calculated  exactly  in  the  high  z  and  low  z 
power  series  expansions  of  6,  the  fraction  of  adsorption  sites  which  are  occupied  by 
adsorbate  (13),  (14).  Using  these  series  expansions,  we  constructed  an  approximant  of  the 
form 


_p(^ 

I  +P(y->) 

which  is  accurate  at  low  z,  and  an  approximan» 


[41 


1 

^  \+P(y) 


[51 


which  is  accurate  at  high  z.  Here  P(y)  is  a  polynomial  in  the  parameter 


[6] 


where  q  is  the  number  of  first  neighbors  to  a  lattice  site  (6  for  the  triangular  lattice  and  3 
for  the  honeycomb  lattice).  The  attractive  lattice  gas,  for  which  w(r)  <  0,  undergoes  a 
first-order  condensation  transition  aty  =  1. 

The  approximants  0^  and  become  equivalent  to  the  Langmuir  adsorption  isotherm 
if  w{r)  =  0,  i.e.,  if  the  lateral  interactions  are  not  present.  These  approximants  are  thus  a 
natural  extension  to  Langmuir's  isotherm. 

Finally,  using  the  switching  function 


n(y.s)  =  ^[\^erf[(y-\)s]]  ,  [7] 

we  construct  a  smooth,  continuous  approximate  isotherm  of  the  form 


0=  r?O’.5)0/+  [1  -  T1(y,s)]df, 


[8] 


which  is  accurate  at  both  high  and  low  z.  The  switching  function  is  used  to  approximate 
the  effects  of  surface  imperfections  and  diffusion  on  the  isotherm.  Certainly,  many 
different  functional  forms  can  be  used  for  the  switching  function  (1 1),  (12). 

A  model  voltammogram  can  be  constructed  from  the  approximate  isotherm  of  Eq.  [8], 
where  a  measure  of  the  current  intensity  is  given  as 

=  191 

dt  d'P  dt 

The  current  intensity  given  by  Eq.  19]  differs  in  sign  from  the  intensity  defined  earlier  (11), 
(12).  Since  the  model  is  completely  reversible,  the  shape  of  the  overall  voltammogram  will 
be  independent  of  the  sign  of  the  intensity,  but  using  /  as  defined  in  Eq.  [9]  is  consistent 
with  the  normal  convention  in  which  the  current  is  positive  during  the  stripping  process, 
during  which  0  is  decreasing. 

A  spike  occurs  in  a  plot  of  I  versus  "Paty  =  1,  the  point  where  the  lattice  gas 
undergoes  a  first-order  phase  transition.  From  Eqs.  [4]-[9]  we  obtain  the  peak  height  at  the 
transition  as 


30; 


[10] 


The  first  term,  which  contains  the  factor  0/,  is  large  for  isotherms  at  temperatures 
well  below  the  critical  temperature  and  vanishes  at  the  critical  point.  The  second  term 
contains  the  factor  which  is  very  small  at  the  transition  at  temperatures  well 


below  the  critical  temperature  and  is  infinite  at  the  critical  point.  Hence,  for  temperatures 
far  from  the  critical  temperature  of  the  transition,  the  peaJc  height  is  proportional  to  the 
change  in  surface  coverage  at  the  transition.  The  factor  s,  which  is  a  measure  of  the 
sharpness  of  the  switching  function,  is  an  approximate  measure  of  the  degree  of  perfection 
of  the  electrode  surface. 

It  is  interesting  to  note  that  the  current  intensity  also  satisfies  the  relationship 


^  dz  dt 


[11] 


where  Kj  is  the  isothermal  compressibility  of  the  equivalent  two-dimensional  lattice  gas. 
Assuming  a  constant  sweep  rate,  Eq.  [9]  can  be  integrated  to  yield 


Thus  the  area  under  a  spike  on  the  voltammogram  is  proportional  to  the  change  in  surface 
coverage  which  occurs  at  the  transition. 

For  temperatures  well  below  the  critical  temperature,  for  which  the  peak  height  given 
by  Eq.  [lOJ  is  roughly  proportional  to  s,  the  width  of  the  peak  must  thus  be  inversely 
proportional  to  s,  because  the  area  under  the  peak,  given  by  Eq.  [12],  is  independent  of  s. 
Hence,  a  rather  broad  peak  will  result  in  the  model  voltammogram  for  systems  in  which 
diffusion  effects  or  surface  imperfections  affect  the  sharpness  of  the  phase  transition. 

Equation  [1]  at  25®C  yields 


r  =  38.922  V  , 


[13] 


where  V  is  the  potential  in  volts. 

’  we  see  from  Ei 
equal  to  I  ^Vj)  -  | . 


In  a  voltammogram  plot  of  /,  in  units  of  I  dVIdt  I 


versus  V'  we  see  from  Eq.  [12]  that  the  area  of  the  voltammogram  from  Vj 


to  is 


AN  EXAMPLE  VOLTAMMCXJRAM 

In  the  presence  of  sulfate  ions,  the  underpotential  deposition  of  copper  on  the  (1 1 1) 
surface  of  gold  produces  a  voltammogram  (3)  which  contains  two  sharp  spikes  with  areas 
in  the  approximate  ratio  of  2: 1 .  This  can  be  interpreted  as  resulting  from  a  sequence  of  two 
first  order  phase  transitions,  corresponding  to  copper  coverages  of  2/3  and  1  monolayer, 
respectively. 


We  have  proposed  (12)  that  this  voltammogram  can  result  due  to  the  following 
sequence  of  events.  First,  the  sulfate  ions  occupy  the  sites  of  the  triangular  lattice  A 
corresponding  to  the  adsorption  sites  on  the  (111)  surface  of  gold.  Because  of  packing 
considerations,  a  pair  of  first-neighbor  sites  cannot  be  simultaneously  occupied  by  sulfate 
ions.  As  a  result,  the  system  undergoes  an  order-disorder  transition,  equivalent  to  that 
which  occurs  in  the  hard  hexagon  lattice  gas  (15),  in  which  one  of  the  three  triangular 
sublattices  Aj  of  A  is  preferentially  occupied  by  sulfate  ions.  This  transition,  which  is 
second  order,  would  produce  a  spike  with  an  infinitesimal  area  which  would  thus  not  be 
observable  as  a  spike  in  the  experimental  voltammogram. 

As  the  potential  is  made  less  positive,  and  nearly  all  the  sites  of  Aj  are  occupied  by 
sulfate  ions,  the  system  will  undergo  a  first-order  phase  transition  as  the  copper  ions  are 
adsorbed  onto  the  vacant  sites  of  A,  which  form  a  honeycomb  lattice  Ah-  Following  the 
development  of  the  previous  section,  the  contribution  to  the  adsorption  isotherm  due  to  this 
process  is  given  as 


where  </  =  3  in  Eq.  [6],  is  the  stickiness  parameter  at  the  sites  of  Au,  and  the  lateral 
interaction  energy  is  Md),  d  being  the  first  neighbor  separation  in  Aj^  anoin  A. 

As  the  potential  becomes  less  positive,  the  copper  ions  then  replace  the  sulfate  ions  and 
occupy  the  sites  of  the  triangular  lattice  Aj,  thereby  completing  the  formation  of  a  copper 
monolayer  on  A.  The  portion  of  the  adsorption  isotherm  due  to  this  process  is  given  by 
our  model  as 


9j  =  tj(y-p,  ix)%'  /  +  [1  - 


where  q  =  6in  Eq.  [6],  Xj  is  the  stickiness  parameter  at  the  sites  of  Aj,  and  w(VJ’d)  is 
the  interaction  energy  between  copper  atoms  on  neighboring  sites  of  Ap 

The  full  adsorption  isotherm  is  then  given  as 

If  the  temperature  is  well  below  the  critical  temperatures  for  the  first  order  transitions  on  A^ 
and  Aj,  then  the  resulting  voltammogram  will  contain  two  spikes  with  areas  in  the  ratio  of 
2:1.  The  resulting  model  voltammogram,  a  plot  of  /  versus  V  where  /  is  in  units  of 
I  dV/dt  I ,  is  illustrated  in  Fig.  1  for  the  choice  of  parameters 


^  =  2. 

ArpO(0, 0)  =  9.1  , 


AhP®(0,  0)  =  5.2  X  105 . 
g^id)  =  14 . 
g^{iTd)  =  3A  . 
j-p  =  5  , 

=  2  . 

This  model  voltammogram  (12)  is  quite  similar  to  the  one  which  has  been  obtained 
experimentally  (3). 
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